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Abstract 

We consider the MSSM with see-saw mechanism of neutrino mass generation and 
soft SUSY breaking with flavour-universal boundary conditions at the GUT scale, in 
which the lepton flavour violating (LFV) decays u — > e + 7, r — > fi + 7, etc., are 
predicted with rates that can be within the reach of present and planned experiments. 
These predictions depend critically on the matrix of neutrino Yukawa couplings Y v 
which can be expressed in terms of the light and heavy right-handed (RH) neutrino 
masses, neutrino mixing matrix C/pmnSj an d an orthogonal matrix R. We investigate 
the effects of Majorana CP-violation phases in C/pmns> an d of the RG running of light 
neutrino masses and mixing angles from Mz to the RH Majorana neutrino mass scale 
Mr, on the predictions for the rates of LFV decays fj, — > e+7, r — > /U+7 and r — > e+7. 
Results for neutrino mass spectrum with normal hierarchy, values of the lightest i/-mass 
in the range < m\ < 0.30 eV, and quasi-degenerate heavy RH Majorana neutrinos in 
the cases of R = 1 and complex matrix R are presented. We find that the effects of the 
Majorana CP-violation phases and of the RG evolution of neutrino mixing parameters 
can change by few orders of magnitude the predicted rates of the LFV decays /i — ► e+7 
and r — > e + 7. The impact of these effects on the r — > u + 7 decay rate is typically 
smaller and only possible for mi ^ 0.10 eV. If the RG running effects are negligible, in 
a large region of soft SUSY breaking parameter space the ratio of the branching ratios 
of the u — ► e + 7 and r — > e + 7 (r — * fx + 7) decays is entirely determined in the case 
of R = 1 by the values of the neutrino mixing parameters at Mz ■ 
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1 Introduction 



The experiments with solar, atmospheric, reactor and accelerator neutrinos [1-5] have 
provided during the last several years compelling evidence for the existence of neutrino 
oscillations caused by nonzero neutrino masses and neutrino mixing. The //-oscillation data 
obtained in the solar neutrino, atmospheric neutrino, KamLAND (KL) and K2K experiments 
imply the existence of 3-u mixing in the weak charged- lepton current (see, e.g., [6]): 



where un are the flavour neutrino fields, is the field of neutrino Vj having a mass rrij and 
U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix [7], U = £/pmns- It 
follows from the results of the neutrino oscillation and 3 H /3-decay experiments [8] that the 
neutrino masses rrij are much smaller than the masses of the charged leptons and quarks: 
rrij < 2.3 eV (95% C.L.) \ 

A natural explanation of the smallness of neutrino masses is provided by the see-saw 
mechanism of neutrino mass generation [11]. The see-saw mechanism predicts the light 
massive neutrinos Vj to be Majorana particles. An integral part of the mechanism are the 
heavy right-handed (RH) Majorana neutrinos. Being singlets with respect to the Standard 
Theory gauge symmetry group [12], the RH neutrinos can acquire Majorana masses which 
are much larger than the electroweak symmetry breaking scale. The CP-violating decays 
of the heavy RH Majorana neutrinos in the Early Universe could generate, through the 
leptogenesis scenario, the observed baryon asymmetry of the Universe [13]. 

In grand unified theories (GUT) the masses of the heavy RH Majorana neutrinos are 
typically by a few to several orders of magnitude smaller than the scale of unification of the 
electroweak and strong interactions, Mqut ~ 2 x 10 16 GeV. The presence of the electroweak 
symmetry breaking scale and of the scale of RH neutrino masses Mr in a theory 2 can lead to 
a hierarchy problem. In SUSY (GUT) theories this problem could be avoided: the hierarchy 
between the two mass (energy) scales is stabilised. Thus, the SUSY (GUT) theories with see- 
saw mechanism provide an internally consistent framework for generation of light neutrino 
masses rrij and of the baryon asymmetry of the Universe. 

The existence of the flavour neutrino mixing, eq. implies that the individual lepton 
charges, L h I = e, /x, r, are not conserved (see, e.g., [14]). Thus, lepton charge non- conserving 
processes like \i~ — > e~ + 7, f/T — > e~ + e + + e _ , r~ — > e~ + 7, — > fi~ + 7, il~ + (A, Z) — > 
e~ + (A, Z), etc. are allowed and should take place. However, if the neutrino mixing in the 
weak charged lepton current, eq. ((H), is the only source of L e , and L T non-conservation, 
as in the minimally extended Standard Theory with massive neutrinos [15], the rates and 
cross sections of the indicated lepton flavour violating (LFV) processes are suppressed by the 
factor [15] (see also [16]) {rrij/M w ) A < 6.7 x 10~ 43 , M w being the W ± mass. This renders 

^^More stringent upper limit on rrij follows from the constraints on the sum of neutrino masses obtained 
from cosmological/astrophysical observations, namely, the CMB data of the WMAP experiment combined 
with data from large scale structure surveys (2dFGRS, SDSS) [9]: J2 3 m j < (0.7-2.0) cV (95% C.L.), where 
we have included a conservative estimate of the uncertainty in the upper limit (see, e.g., [10]). 

2 In the case of three heavy RH Majorana neutrinos Nj with masses Mj we will consider, Mr is determined 
by the mass of the lightest RH Majorana neutrino, Mr = min(Mj). 
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them unobservable. Indeed, the existing upper limits on, e.g. /x — > e + 7, /x — > 3e, and 
r — > /x + 7 decay branching ratios and on the relative cross section of the /x — e conversion 
process /x~ + Ti — > e~ +Ti, are at least by 30 orders of magnitude bigger [17-19] (90% C.L.): 

BR(/x -»■ e + 7) < 1.2 x 10" 11 , BR(zx -> 3e) < 1.2 x 10~ 12 , 

(2) 

BR(r -> /x + 7) < 6.8 x 10" 8 , R(/x~ + Ti -> e~ + Ti) < 4.3 x 10~ 12 . 

Future experiments with increased sensitivity can reduce the current bounds on BR(/i — > 
e + 7), BR(r — > /x + 7) and on R(/x~ + (A, Z) — »■ e~ + (A, Z)) by a few orders of magnitude 
(see, e.g., [20]). In the experiment MEG under preparation at PSI [21] it is planned to reach 
a sensitivity to 

BR(/x-^e + 7) ~ (10~ 13 - 10" 14 ) . (3) 

It has been noticed in 1986 that in SUSY (GUT) theories with see-saw mechanism of 
neutrino mass generation, the rates and cross sections of the LFV processes can be strongly 
enhanced [22] and could be close to existing upper limits. If the SUSY breaking occurs 
via soft terms with universal boundary conditions at a scale Mx above the RH Majorana 
neutrino mass scale 3 Mr, Mx > Mr, the renormalisation group effects transmit the LFV 
from the neutrino mixing at Mx to the effective mass terms of the scalar leptons at Mr, 
generating new LFV corrections to the flavour- diagonal mass terms. For scalar lepton masses 
of ~ few hundred GeV, the LFV mass corrections at Mr are typically of the order of a few 
GeV and thus are much larger than the light neutrino masses rrij. As a consequence (and 
in contrast to the non-supersymmetric case), the LFV scalar lepton mixing at Mr generates 
additional contributions to the amplitudes of the LFV decays and reactions which are not 
suppressed by the small values of neutrino masses rrij. As a result, the LFV processes can 
proceed with rates and cross sections which are within the sensitivity of presently operating 
and future planned experiments [22,24] (see also, e.g., [25-31]). 

One of the basic ingredients of the see-saw mechanism is the matrix of neutrino Yukawa 
couplings, Y u . Leptogenesis depends on Y„ as well. In the large class of SUSY models with 
see-saw mechanism of neutrino mass generation and flavour-universal soft SUSY breaking at 
a scale Mx > Mr, we will consider, the probabilities of LFV processes also depend strongly 
on Y v (see, e.g., [25,26]). Thus, the matrix of neutrino Yukawa couplings Y„ connects in 
the see-saw theories the light neutrino mass generation with leptogenesis; in SUSY theories 
with soft flavour-universal SUSY breaking at Mx > Mr, Y u links the light neutrino mass 
generation and leptogenesis with LFV processes (see, e.g., [28,29]). 

The matrix Y„ can be expressed in terms of the light neutrino and heavy RH neutrino 
masses, the neutrino mixing matrix ?7pmns, and an orthogonal matrix R [25]. Leptogenesis 
can take place only if R is complex. Obviously, Y " v depends on the Majorana phases in the 
PMNS matrix £/"pmns- It was shown in [28] (see also [26]) that in the case of R ^ 1, the 
Majorana CP-violation phases in C/pmns [32] can affect significantly the predictions for the 
rates of LFV decays /x — > e + 7, r— > /x + 7, etc. in the class of supersymmetric theories of 
interest . 

The matrix ~Y V can be defined, strictly speaking, only at scales not smaller than Mr. 
The probabilities of LFV processes depend on Y„ at the scale Mr, Y v = Y„(Mr). In order 

3 The possibility of "flavour-blind" SUSY breaking of interest is realised, e.g. in gravity-mediated SUSY 
breaking scenarios (see, e.g., [23]). 
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to evaluate Y u (Mr) one has to know, in particular, the light neutrino masses rrij and the 
mixing matrix £/pmns a ^ Mr, i.e., one has to take into account the renormalisation group 
(RG) "running" of rrij and C/pmns from the scale M z ~ 100 GeV, at which the neutrino 
mixing parameters are measured, to the scale Mr [33-37]. 

Working in the framework of the class of SUSY theories with see-saw mechanism and 
flavour-universal soft SUSY breaking at a scale Mx > Mr, we investigate in the present 
article the effects of RG running of neutrino mixing parameters and of the Majorana CP- 
violation phases on the predictions for the rates charged lepton flavour violating decays. We 
consider neutrino mass spectrum with normal hierarchy (NH), m 1 < m 2 < m 3 , allowing 
the lightest neutrino mass mi to vary in the range < mi < 0.30 eV. The latter covers 
all possible varieties of light neutrino mass spectrum of the NH type: normal hierarchical, 
quasi- degenerate as well as spectrum with partial normal hierarchy (see, e.g., [6]). In order 
to find the changes of the neutrino masses rrij and mixing angles in C/pmns due to the RG 
evolution from the scale M z ~ 100 GeV to M R , the corresponding 1-loop RG equations 
for the see-saw Majorana mass matrix of the left-handed neutrinos were solved numerically. 
The best fit values of the known solar and atmospheric neutrino mixing angles and neutrino 
mass squared differences at Mz were used as "initial conditions" for the solution of the RG 
equations. The values of the neutrino mixing parameters at the scale Mr served as input in 
the calculation of Y„(Mr) and correspondingly of the rates of the LFV decays Z,; — > lj + 7. 
The latter were performed for quasi-degenerate heavy RH Majorana neutrinos for a number 
of typical sets of values of the soft SUSY breaking parameters. We have obtained results 
both in the cases of R = 1 and of complex matrix R. 

Detailed predictions for the rates of the LFV processes in the class of SUSY models 
with see-saw mechanism considered were obtained, e.g. in refs. [25-31]. However, our study 
overlaps little with those performed in these articles. 



2 Neutrino Mixing Parameters from Neutrino Oscillation Data 



We will use the standard parametrisation of the PMNS matrix C/pmns ( see > e -9-i [38]): 
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is the Dirac CP- 



violating phase and a and /3m are two Majorana CP- violation phases [32,39]. One can 
identify the neutrino mass squared difference responsible for solar neutrino oscillations, Am? 
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> 0. The neutrino mass squared difference driving 



oscillations of atmospheric i/„ 



is then given by 



^> Am^. The corresponding solar and atmospheric neutrino 
mixing angles, 6 & and 9 a, coincide with 6*12 and 6*23, respectively. The angle #13 is limited 
by the data from the CHOOZ and Palo Verde experiments [40]. 

The existing neutrino oscillation data allow us to determine Am| 1; jAm^l, sin 2 9 12 and 
sin 2 2#23 with a relatively good precision and to obtain rather stringent limits on sin 2 #13 
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(see, e.g., [41,42]). The best fit values and the 95% C.L. allowed ranges of Am^, sin 2 #i 2 , 
[Am 2 ^ and sin 2 2# 23 read 4 : 

Am^ = 8.0 x 10~ 5 eV 2 , sin 2 21 = 0.31 , 

(5) 

Am^ = (7.3 - 8.5) x 10~ 5 eV 2 , sin 2 9 12 = (0.26 - 0.36) , 
I Am 2 I = 2.2 x 10~ 3 eV 2 , sin 2 29 23 = 1.0 , 

I 311 

lAmfJ = (1.7-2.9) x 10~ 3 eV 2 , sin 2 2# 23 > 0.90 . 
A combined 3-v oscillation analysis of the solar neutrino, KL and CHOOZ data gives [41,42] 

sin 2 fl 13 < 0.027 (0.047), at 95% (99.73%) C.L. (7) 

The neutrino oscillation parameters discussed above can (and very likely will) be measured 
with much higher accuracy in the future (see, e.g., [6]). 

The sign of Am 2 A = Am^, as it is well known, cannot be determined from the present (SK 
atmospheric neutrino and K2K) data. The two possibilities, Am 2 ^^ > or Am^^j < 
correspond to two different types of z/-mass spectrum: 

- with normal hierarchy m\ < m 2 < m 3 , Am\ = Am| x > 0, and 

- with inverted hierarchy m 3 < mi < m 2 , Am\ = Am| 2 < 0. 

Depending on the sign of Am\, sgn(Am 2 v ), and the value of the lightest neutrino mass, 
min(mj), the i^-mass spectrum can be 

- Normal Hierarchical: m-i-C m 2 <Cm 3 , m2=(Amg)^~ 0.009 eV, m 3 =|Am^|5~ 0.045 eV; 

- Inverted Hierarchical: m 3 <C mi < m 2 , with m,i j2 = |Am||2 ~ 0.045 eV; 

- Quasi- Degenerate (QD): m\ = m 2 = m 3 = tuq, m 2 ^> jAm^], niQ <; 0.10 eV. 

It is well-known, the theories employing the see-saw mechanism of neutrino mass gener- 
ation [11] we are going to discuss further, predict the massive neutrinos Uj to be Majorana 
particles. Determining the nature of massive neutrinos is one of the most formidable and 
pressing problems in today's neutrino physics. Being of fundamental importance for under- 
standing the origin of ^-masses and mixing, it is widely recognised as a highest priority goal 
of the future studies of neutrino mixing (see, e.g., [6,45]). 

If the massive neutrinos Vj are Majorana fermions, getting information about the Majo- 
rana CP-violation phases in ?7pmns would be highly problematic. The oscillations of flavour 
neutrinos, v\ —>■ vy and v\ — > vy , 1,1' = e,fi,T, are insensitive to the Majorana CP- violation 
phases a and (5m [32,46]; they are insensitive to the nature - Dirac or Majorana, of massive 
neutrinos Uj. The only feasible experiments that at present have the potential of establish- 
ing the Majorana nature of light neutrinos z/j and of providing information on the Majorana 
CP- violation phases in ?7pmns are the experiments searching for the neutrinoless double beta 
(((3(3)ov-) decay, (A, Z) -> (A, Z + 2) + e~ + e~ (see, e.g., [14,45,47]). If the (/^V-decay is 
observed, the measurement of the (/3/3)oi/-decay half-life combined with information on the 
absolute scale of neutrino masses (or on min(mj)), might allow to significantly constrain, or 
possibly determine, the Majorana phase a [38,48], for instance. Establishing CP- violation 
due to the Majorana phases in ?7pmns, however, was found to be remarkably challenging [49]. 

4 The data imply, in particular, that the high-LMA solution (see, e.g., [43]) is excluded at ~ 3.3ct. Maximal 
solar neutrino mixing is ruled out at ~ 6a; at 95% C.L. one finds cos2#q > 0.28 [41], which has important 
implications [44]. 
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3 The See-Saw Mechanism and Neutrino Yukawa Couplings 



We consider the minimal supersymmetric standard model with RH neutrinos (MSSMRN). 
The superpotential for the lepton sector is given by 

m =£?(Y E ) ii L j • H d + NftYJijLj ■ H u - J(M N )^JV| , (8) 

where Lj, and iVJ, j = e, [M, t = 1,2,3, are the chiral supermultiplets respectively of 
the SU(2) L doublet lepton field Lj L and of the SU(2) L singlet charged lepton and neutrino 
fields lj L = C(ljn) T and N? L = CNj R) where Ijr and Nj R are RH fields and C is the charge 

conjugation matrix, while H u and Hd are the supermultiplets of the two Higgs doublet fields 
H u and Hd carrying weak hypercharges — | and |. In eq. (JHJ), Ye is the 3x3 matrix of 
the Yukawa couplings of the charged leptons, Y v is the 3x3 matrix of neutrino Yukawa 
couplings, and M« is the Majorana mass matrix of the RH neutrinos Njr. 

In the framework of MSSMRN one can always choose a basis in which both Ye and Mn 
are diagonal. Henceforth, we will work in that basis and will denote by D« the corresponding 
diagonal RH neutrino mass matrix, Dn = diag(M 1; M 2 , M 3 ). 

Below the see-saw scale, Mr = min(Mj), the singlet supermultiplets containing the 
RH neutrino fields are integrated out, and the superpotential, eq. (JEJ), becomes 

Wi =Kf(Y E )^4- • H d - -(KnY^U ■ H u )(Lj ■ H u ) , (9) 

where 

= {YjY\M^) kl {Y u f . (10) 

The last term yields Majorana mass term for the left-handed (LH) flavour neutrinos after 
the electroweak symmetry breaking, 

C v m = -\v$ j {m v y h v L k + h.c., (11) 

where v^- = C{pLj) T and 

(m v ) ij = {K N f{H u f . (12) 

Here (H u ) is the vacuum expectation value of H u , (H u ) =v u = v sin ft, where v = 174 GeV. 
The neutrino mass matrix is related to the light neutrino masses rrij and the PMNS mixing 
matrix as 

{m v f = (U*) ik m k (U^ . (13) 

Using ()12j1 and (JEJ), we can rewrite the matching condition, eq. (JTUJl . at the energy scale 
Mr in the form 

U*(Mr) D„(M r ) V\M r ) = -v 2 u YI(Mr)M^(M r )Y u (Mr) . (14) 
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where D„ = diag(m 1; m 2 , m 3 ). Thus, in the basis in which the RH neutrino mass matrix is 
diagonal M N = D N , the neutrino Yukawa coupling at M R can be parametrised as [25] 

Y V (M R ) = —y/B N (M R ) R y/B u (M R ) TJ*(M R ) . (15) 

Here R is a complex orthogonal matrix, R T R = 1. The PMNS matrix U and the light 
neutrino masses in D„ in eq. (j!5|) should be evaluated at M R with the help of the renormal- 
isation group equations which describe their evolution from the the scale Mz, where they 
are measured 5 , to the scale M R . 

4 The Renormalisation Group Evolution Effects 

The RG running corrections to the neutrino mass matrix in MSSM have been studied 
by many authors (see, e.g., [36,37] and the articles quoted therein). In order to find the 
RG evolution from Mz to M R of the neutrino mixing angles 9u, #23, #13 and the masses rrij, 
we have solved numerically the relevant 1-loop RG equation for the neutrino mass matrix, 
eq. (fT2"|). in MSSMRN [33] (see also, e.g., [30,36,37]). The case of neutrino mass spectrum 
with normal hierarchy at Mz, mi(Mz) < m 2 (Mz) < m 3 (Mz), was considered. The best 
fit values of the neutrino oscillation parameters #12, O23, Am^ and Am^ at Mz, given in 
eqs. (jSJ and (JBJ), were used as input in the calculation of the RG running effects. The CHOOZ 
mixing angle at Mz was set to zero, sin #13 (Mz) = 0. 

We will summarise below some of the results regarding the RG running of the neutrino 
mixing parameters, which are relevant for our further analysis. 

• The corrections can become significant for relatively large tan /5. 

• The mixing angles are stable with respect to the RG running in the case of normal 
hierarchical neutrino mass spectrum, mi <C m 2 <C m 3 . 

• For mi ^ 0.05 eV, the RG running effects can change significantly the neutrino (lepton) 
mixing angles provided tan /3 is sufficiently large. 

• The solar neutrino mixing angle 6 Q = #i 2 at M R depends strongly on the Majorana 
phase a, which plays very important role in the predictions of the effective Majorana 
mass in (f3(3) 0u — decay (see, e.g. [47]). The effect of RG running for 12 is largest when 
a takes the CP-conserving value a = 0; it is relatively small or negligible for the 
alternative CP-conserving values of a [50], a = ±7r. For a = and tan/5 ~ 50, for 
instance, we have tan 2 ^i 2 (M R ) < 0.5 tan 2 Qyi(M z ) for mi > 0.02 eV. 

• The RG running effects depend weakly on the Majorana phase (3m- 

• The atmospheric neutrino mixing angle &a = #23 is essentially stable against the RG 
running corrections, except in the cases when mi and/or tan/3 are/is relatively large, 
i.e., (i) when mi ^ 0.2 eV and a ^ 37r/4 if tan/5 ~ 10, and (ii) for any m 1 and a if 
tan/5 > 40. 

5 Strictly speaking, the neutrino mixing parameters are measured at scales ~ 1 GeV and ~ 1 MeV which 
are smaller than Mz- However, these physical parameters do not get any significant corrections from the 
RG running below Mz- 



6 



• The correction to sin # 13 is relatively small even for large tan (3 if mi ^ 0.05 eV, and for 
any mi ^ 0.30 eV if 9i 3 (M z ) = and a = 0. For a = ir and tan/3 ~ 50 one can have 
sm9 13 (M R ) > 0.10 formx > 0.08 eV even if sin 9 13 (M Z ) = 0. If, e.g. sin# 13 (M z ) = 0.10 
eV, the RG effects lead for a = (tt/2 ^ a ^ ir) to considerably smaller (larger) values 
of sin6»i 3 (M J? ) at m 1 > 0.05 eV; for a ~ tt/4 one has sin6 13 (M z ) = sin9 13 (M R ). 

• For sufficiently large tan/3, tan/3 ^ 30, the value of Am| 1 (M fi ) depends strongly on 
mi m the interval rrii ^ 0.05 eV, and on a for mi ^ 0.10 eV. In the case of tan/3 ~ 
50, for instance, we have Am|(M fl ) = 3Am 2 21 (M z ) (Am 2 21 (M R ) = 10Am 2 21 (M z )) at 
mi = 0.05 (0.10) eV. The dependence of Am 31 (M#) on mi and a is rather weak, unless 
tan (3 > 40 and mi > 0.10 eV. 

• Some products of the neutrino mixing parameters, such as Si2Ci2C2 3 {mi — m 2 e ta ) are 
practically stable with respect to RG running (for further discussion see Appendix A 
and, e.g. refs. [30,36,37]). 

Similar results hold also in the case of neutrino mass spectrum with inverted hierarchy. 

The RG equation for the neutrino mass matrix in MSSMRN and its approximate solution 
as well as analytic expressions for some of the neutrino mixing parameters at M R in the case 
of interest are given in Appendix A. The specific features of the effects of RG running of 
the neutrino mixing angles and mass (squared) differences discussed above are illustrated in 
Figs. [U and El where the values of tan 2 9i 2 , tan 2 # 2 3 and sin^, and of (m 2 — mi), Am^, 
(m 3 — mi) and Am 31 , at M R = 2 x 10 13 GeV are shown as functions of mm(mj) = mi. The 
neutrino mass spectrum is assumed to be with normal hierarchy at M z . The two figures 
correspond to the following values of the neutrino oscillation parameters and (3m at M z : 
tan 2 612 = 0.40, tan 2 9 23 = 1, Am 2 21 = 8.0 x 10~ 5 eV 2 , Am 2 31 = 2.2 x 10~ 3 eV 2 , sin#i 3 = 
and Pm = 0. 

The results obtained in the present Section are used in the calculations of the U — > lj + 7 
decay rates to which we turn next. 

5 Lepton Flavour Violation in MSSMRN and k — > lj ; + 7 Decays 

In the MSSMRN under discussion, the boundary conditions for the soft SUSY breaking 
are assumed to take place at some scale Mx > Mr, typically at Mx ~ Mqut- They are 
described in the general case by the Lagrangian 

Aoft = - l -M g3 G^ c G^ - ^MgzWfWi - \m 9 iB c B 

- (M 2 L y%l Lj - {M 2 E f~e R fe R] - {{A E f~e*J Lj ■ H d + h.c) 

- (m 2 m \H d \ 2 - m 2 H2 \H u \ 2 ) - (BfiH d ■ H u + h.c.) (16) 

where G (a) , W u and B are the gauginos associated with the SU(3) C , SU(2)^, and U(l)y 
gauge symmetry groups, and H u , H d , l Li , and e* Ri are the scalar components of H u , H d , Lj, 
and E? superfields, respectively, and the dots mean additional soft RH sneutrino and squark 
terms which do not play a role in our further discussion. We will consider in what follows the 
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class of MSSMRN in which the soft SUSY breaking terms in eq. (|16|) are flavour-universal 
at the scale M x , 



M gi = mi/2 , 

m Hl = m H2 = m l 7 

{M L f = {M E f = ml 8 ij , 

{A E f = A {Y E f, A = a m , (Y E f = Y E ^. (17) 

In this case the only lepton flavour mixing present at Mx is that due to the neutrino Yukawa 
couplings. However, at the scale M R = min(Mj) = M 1; where the RH neutrinos decouple, 
the RG running effects induce flavour non-diagonal elements in the effective soft SUSY 
breaking Lagrangian. The off-diagonal elements of interest for our further discussion are 
given in the leading- log approximation by [24], 

(MIT * - i™o(3 + kl 2 ) In ^ • + j, (18) 



S7T 



k 



{M%)V~0, i^j, (19) 



{A E ) ij * - ^rn a Y E ^ ¥f Y» j , Mi- 

k 



These terms generate new contributions in the amplitudes of the LFV processes. In the 
"mass insertion" (see, e.g., [24,26,30]) and leading-log approximations, the branching ratio 
of li — > lj + 7 decay due to the new contributions has the following form 



T(/j — > evv) a 



3 



cm 



total (h) Gpm? s 



(3 + al)ml 



In 2 



^ K v) M k v 

k 



tan 2 /?, (21) 



where mj represents SUSY particle mass. It was shown in [30] that in most of the relevant 
soft SUSY breaking parameter space, the expression 

m% ~ 0.5 ml m\ /2 (m 2 + 0.6 m\ /2 f , (22) 

gives an excellent approximation to the results obtained in a full renormalisation group 
analysis, i.e., without using the leading- log and the mass insertion approximations. Thus, 
in the case of soft SUSY breaking mediation by flavour- universal terms at Mx > Mr, the 
predicted rates of LFV processes such as \i — > e+7 decay are very sensitive to the off-diagonal 
elements of 

Yt(M R )Y v (M R ) = 1 U(M R )x/n u (M R ) Rt B N R \/~D u (M R )\j\M R ) . (23) 

V u 

Obviously, Y^Y^ is affected by the RG running of U and D„ below the scale M R , which is 
strongly dependent on tan/5, the light neutrino mass m\ and the Majorana phase a. 
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5.1 Results for R = 1 



In our further analysis we consider the case of quasi-degenerate heavy RH neutrinos, 
Mx = M 2 = M 3 = M R and we set M x = M GUT = 2 x 10 16 GeV. We also first fix R in 
eq. ()15j) to be the unit matrix, R = 1. The case of R ^ 1 is discussed later. Under the 
conditions specified above the branching ratios BR(/j — > lj + 7) depend on Yj Y^, which in 
turn takes the form 

Mr 

Yt(M R )Y u (M R ) = -^V(M R )dia,g(m 1 (M R ),m 2 (M R ),m 3 (M R ))V\M R ) . (24) 

v 

u 

For the off-diagonal elements of Yj Y v of interest we get to leading order in small parameters 



(K^)l2 — ^2lC23 c 12-5l2Cl3 + A 31 •523 c 13 s 13 e 

16 + 0(A 21 a 13 ) , (25a) 

(Yjr^g = - A 21 s 2 3C 12 s 12 c 13 + A 31 c 23 c 13 s 13 e- tS + 0(A 21 s 13 ) , (25b) 

(Y}Y U ) 23 = A 31 s 23 c 23 c 2 13 + 0(A 21 , A jlSl3 , s 2 13 ) , (25c) 

where all angles should be evaluated at M R and 

Mr Mr Am? A Mr) 

Aij = -4 (mAM R ) - mj (M R )) = lJ 2 (m \ ■ ( 26 ) 

vl vl mi{M R ) + mj{M R ) 

5.1.1 The Case of Negligible RG Corrections 

The expressions for (YJY^m, k ^ I = 1, 2, 3, we have derived, eqs. (|25aJ) - (|25cj) . are 
valid both for light neutrino mass spectrum with normal and inverted hierarchy, and, corre- 
spondingly, for NH, IH and QD spectrum. If in the case of QD heavy RH neutrinos the RG 
effects on the angles % and masses mj are negligible or sufficiently small (e.g. if tan/3 ^ 10), 
in the large region of the soft SUSY breaking parameter space where the approximate ex- 
pression ()21|) with m$ given by eq. (J2*2*|) reproduces with a relatively high precision the values 
of BR(/j — > lj + 7) calculated exactly, the ratios 

R(ij/km) = ™ ( /' - ) + 7) . , i^j,k^m, (27) 

depend only on the neutrino mixing parameters 6\ 2 , 6 23 , 6i 3 , 5 and mj at Mz- We have, for 
instance, in the case of NH neutrino mass spectrum: 



BR(/i -> e + 7) 



\/ Am 2 21 c 23 ci 2 s 12 + ^ Am 2 31 s 23 s 13 e~ 



R(21/31) = —j- ^ — 2 , NH . (2 



BR(r -> e + 7) 



\/ Am 2 2l s 23 ci 2 si 2 - ^/AmJ^c 23 s 13 e~ 



The expressions for IH (QD) spectrum can obtained by replacing y'Am^ with (— 2y / |A7n||~[) 
(^/Am^ and a/ Am 31 with Am^ and Am 31 ) in eq. (|2*%|) . Depending on whether the terms 
oc \J Amjj, or the terms cx ^/lAm^l, dominate, we get independently of the type of neutrino 
mass spectrum 

R(21/31) ^ tan^ 2 9 23 or tan 2 9 23 . (29) 
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In both cases BR(/i — > e + 7) and BR(r —> e + 7) are predicted to differ at most by a factor 
~ 5, which reflects the current 3<r uncertainty in the value of tan 2 #23- 
Similarly, we have for the ratio R(21/32): 



NH 



y i\m 21 c 12 s 12 s 23 -\- y i\m 31 s 13 c 23 e , in. 

(30) 



The expression for IH (QD) spectrum can be obtained again by replacing y 7 Am^ with 
(— 2y / |Am| 1 |) (y/ Am 2l and a/ Amf^ with Am 21 and Am^) in eq. (1301) . Now the dominance 



of the terms oc y Am^i implies, e.g. for IH spectrum 



R(21/32)= , A A m |, C ¥ 2 , IH - (31) 
4|Am§ 1 | sl 3 



In the case of NH spectrum this ratio is 4 times bigger, while for QD spectrum the fac 
tor (Am2 1 /(4|Am| 1 |)) should be replaced by (Am 2l /\Am 31 \) 2 . If the terms oc y 7 ] A 



311 



dominate, we have independently of the type of the spectrum: 

R(21/32) = s\ 3 c 23 2 . (32) 

In any of these cases we get BR(/i — > e + 7) <C BR(r — > // + 7). 

Let us add that for the current best fit values of the neutrino oscillation parameters 
and sin 2 ^13 close to its current upper limit ((Zj), the terms oc y 7 Am 2l and oc y / ]Am 3 \| in 
R(21/31) and R(21/32) are of the same order of magnitude for NH and IH neutrino mass 
spectra, while the term oc lAm 2 ^ is the dominant one for QD spectrum. Thus, if sin 2 9± 3 
is sufficiently large and the neutrino mass spectrum is of the NH or IH type, the ratios 
R(21/31) and R(21/32) will depend on the Dirac CP-violation phase as well. 



5.1.2 RG Effects 

We turn next to the discussion of the RG effects. In order to find the changes due to 
the RG running from Mz to Mr of the neutrino mass matrix and correspondingly of the 
neutrino mixing angles 9 U , 9 23 , # 13 and masses rrij, the 1-loop RG equations in MSSMRN 
for the matrix in eq.(10), as we have already discussed, are solved numerically. As 
in the previous Section, the neutrino mass spectrum at Mz is assumed to be with normal 
hierarchy, m\{Mz) < TU2(Mz) < m^{Mz)- The heavy RH neutrino mass scale Mr is set 
at Mr = 2 x 10 13 GeV. The same best fit values of the neutrino mixing parameters at Mz 
were used as input in the calculation of the RG running effects. In particular, we have fixed 
sin 9i 3 = at Mz ■ The values of the neutrino mixing parameters at the scale Mr were used 
as input in the calculation of the rates of the LFV decays U — > lj + 7, which were performed 
for values of tan/3 from the interval ~ (10 — 50) and a number of typical sets of values 
of the soft SUSY breaking parameters m , m^ 2 and A in the few 100 GeV range. These 
calculations have been performed without employing the reading-log and the mass insertion 
approximations. 

The results we have obtained using the methods described above show that if R = 1, 
taking into account the effects of the RG running of the neutrino mixing parameters in 
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the case of relatively large tan/? and/or m 1; can change significantly the predictions for 
the LFV U —> lj + 7 decay rates and the corresponding branching ratios BR(/j — > lj + 7) . 
The magnitude of this change depends strongly on the value of the Majorana CP-violation 
phase a. These conclusions are illustrated in Fig. El where the predicted U — ► lj + 7 decay 
branching ratios in the case of R = 1 are shown as functions of the light neutrino mass 
mi for two typical sets of soft SUSY breaking parameters and three values of the Majorana 
CP-violation phase a, a = 0; tt/2; it. 

We find, in particular, that for 0.05 ^ mi ^ 0.30 eV and 30 ^ tan/5 ^ 50, the predicted 
fi — * e + 7 and and r — > e + 7 decay rates, or BR(/i — > e + 7) and BR(r — > e + 7), can 
be enhanced by the effects of the RG running of % and rrij by 1 to 3 orders of magnitude 
if 7r/4 ^ a ^ 7r, while BR(r — > fi + 7) can be enhanced by up to a factor of 10. The 
enhancement of BR(/i — > e + 7) and BR(r — > e + 7) is due to the "regeneration" of the terms 
oc A 31 s 13 i n (Xj^)i2 and (YjY^^ at M R : these terms are zero at M z since we have set 
siz(Mz) = 0. The corresponding enhancement factor in the case of /i — > e + 7 and r — > e + 7 
decays is given approximately by the ratio 

F ^ ( Ar »li m i+ m 2 £is(Mr) \ 2 ^ 
\Am21 mi + m 3 C12S12 / 

where all quantities, and not only S13, should be evaluated at Mr. For a = it, however, 
#12 {Mr) = 12 (M Z ) (see Fig. |TJ) and F is reproduced with a rather good approximation if 
one uses the values of Am^, Am^ and of rrij at M^. For m x = 0.05 (0.10) eV, for instance, 
we have for tan/3 = 50, si 3 (M^) = 0.05 (0.12) (see Fig. |T} and for the enhancement factor 
we find F = 7 (60). The case of r — > e + 7 decay is somewhat more complicated because of 
the possibility of accidental mutual cancellation for a = it of the two terms in eq. (J25bj) at 
Mr 6 . This can lead to a rather strong suppression (by more than an order of magnitude) 
of BR(r — > e + 7) in rather narrow interval of values of m\. The "accidental" suppression 
under discussion can take place only in a limited region of the soft SUSY breaking parameter 
space. 

For a = 0, the combined effect of the RG running of 812, 623 &13 and rrij on BR(/j — > lj+ r y) 
is relatively small even for mi ~ (0.2 — 0.3) eV and tan/? ~ 50. This is a consequence 
of the fact that for a = 0, S13 essentially does not "run" and sis(Mr) = s\^{Mz), the 
effect of RG running of C23 (^23) is relatively small, while the effects of RG running of A 2i 
compensates the effect of running of C12S12 (see also Figs. HJandEJ). In this case BR(/i — > 
e + 7) and BR(r — > e + 7) (BR(r — > /x + 7)) scale as functions of mi approximately as 
~ [AmLjJ^fffli + rn 2 (3))} 2 , where m 2 (3) = (ml + AmL 31 j)a and the masses rrij and 
Am|i 31 can be taken at Mz- 

The effects of the RG running of the neutrino mixing parameters 9 12 , 623, 0\z an d rrij 
are, in general, considerably smaller for tan/3 ~ 10. They are practically negligible for 
BR(r — > fi + 7). If7r/2^a^7r, they can lead to an enhancement of BR(/i — * e + 7) 
approximately by a factor of 2.5 (20) for mi = 0.10 (0.30) eV and are essentially negligible 
if mi < 0.06 eV. In what regards the r — >• e + 7 decay branching ratio, BR(r — > e + 7), 
also in this case it can be strongly suppressed for a = it due to the accidental cancellation 
discussed above. Such a suppression can be realised for relatively large mi ^ 0.15 eV only. 

6 Note that since we set S13 = at Mz, such a cancellation can take place only due to the fact that the 
RG running can generate a non-zero and sufficiently large S13 at Mr. 
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5.2 The Case of R ^ 1 



We shall analyse next briefly the case of nontrivial matrix R: R ^ 1, Successful 
leptogenesis can take place only if R is complex, so we will consider (R)* 7^ R. The matrix 
R of interest, being a complex orthogonal matrix, can be written in the form: R = O e lA , 
where O is a real orthogonal matrix and A is a real antisymmetric matrix, (A) T = —A. In 
the case of QD heavy RH Majorana neutrinos we can set O = 1 without loss of generality. 
Thus, in the subsequent calculations we use R = e tA , parametrising the real antisymmetric 
matrix A as follows 



(34) 



where a, b, c are real. The following representation of e lA proves convenient for our analy- 
sis [28]: 




„ ■ . cosh r — 1 9 sinh r , , 
R = e jA = l A + i — - — A , (35) 



with r = \J a 2 + b 2 + c 2 . The requirement of successful leptogenesis in the case of QD light 
and heavy RH neutrino mass spectra implies [28] that abc 7^ and that 10~ 5 ^ \abc\ <C 1. 
In what follows we will consider the range of values of the parameters a, b, c determined by 
|a|, |6|, \c\ < 0.10. 

The matrix of neutrino Yukawa couplings in the case under discussion has the form 
Y v = ^ e^VDytf. The off-diagonal elements of Y^Y^ of interest are given to leading 
order in the small quantities by 



M R 



7j2 
11 



(Y}Y V ) & = A21C23C12S12C13 + A 3 iS23Ci3Si 3 e 15 

flvW^ (c23(ci 2 e _i t + s \ % j% ) + s lz ci2S 12 s 2 z{e i{ ^ &) - e" i{ ? +<5) )) 

+ Cy/m^s 23 {s^e 1 ^^ + s 13 c 12 e- i(fL ^ J1£+<5) ) + 0(s 2 13 



0(r 2 ) , (36) 



(YjY^is — — A21S23C12S12C13 + A 31 C23Ci3Si3e lS 

a^/m x m2 (-s 23 (c? 2 e" if + s? 2 e*f ) + s 13 c 12 s 12 C2 3 {e l ^- &) - e^^)) 

r / ( i 0M i(^ M H\\ 

+ by/mim 3 ^Ci2C 23 e 2 - S13.S12.S23e 1 2 ') 

+ c^hS ( S i2C 23 e^ + s 13 c 12 S2 3 e- t ^ +5) ) + 0(4,) + 0(r 2 ) , (37) 



Mr 



V 2 
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C^/^i/)23 — ^31 s 23 c 23 c 13 

i^/mrm (-2zci 2 si 2 C23S23 sin ^ + s 13 [c 2 12 (sl 3 e~ l ^~ s) + c 2 23 e % ^~ &) ) 



Mr 



7j2 



b^m x m z \ -si 2 (s 23 e 2 + c 23 e 2 ) + S13C12C23S23 2zsin( — - b) 
+ Cy/m 2 m 3 c 12 {c 23 e 2 + s 23 e 2 ) - S13S12C23.S23 2i sm ( ^ I- 0) 



0(s 



+ 0(r J ) , (3* 



Equations (J3l)j) - (j3"%j) are valid for any of the possible types of light neutrino mass spec- 
trum. The above results imply that in the absence of significant RG effects, the "double" 
ratios R(21/31) = BR(/i -> e + 7)/BR(r -> e + 7) and R(21/32) = BR(/i -> e + 7)/BR(r -> 
/i + 7), depend in the region of validity of eqs. (}2~T]) and ()22|) in the relevant SUSY parameter 
space, not only on 9^, 5 and rrij at Mz, but also on the Majorana CP-violation phases a 
and (3m and on the "leptogenesis parameters" a, b and c. 

The results we find for BR(7j — > lj + 7) taking into account the RG effects in the case 
of complex R given by eq. (|35j) . are illustrated in Figs. E]andEl Both figures correspond to 
(3m = 0, a = 0; 7r/2; it, and the values of 6>i 2 , #23, 6>i 3 , Am 21 and Am^ at Mz used in Figs. IB- 
El In Fig. |U the branching ratios BR(/i — > e + 7), BR(r — > e + 7), and BR(r — > /i + 7) 
are shown as functions of r = V& 2 + fr 2 + c 2 varying in the interval 10~ 4 — 5 x 10~ 2 for 
fixed mi = 0.06 eV and SUSY parameters tan/3 = 10, mi/ 2 = 250 GeV, m = 100 GeV 
and A = —100 GeV. Figure El is similar to Fig. El but shows results for complex R. It 
is obtained assuming that the three constants a, b, c parametrising the matrix R lie in the 
interval —0.1 < a, b, c < 0.1. For comparison, in Fig. El results for R = 1 are also shown 7 . 

In the case of tan/3 ~ 10 and m-i ~ 0.06 eV (Fig. HJ), the RG effects are negligible. For 
a = 0, the branching ratios BR(/i — > e + 7) and BR(r — > e + 7) can be significantly enhanced 
by the effect of complex R [28] for r <; 5 x 10~ 3 : the enhancement factor can be as large as 
~ 100 at r = 5 x 10~ 2 . Similar enhancement can take place for values of a in the interval 
[7r/2,7r]. However, as Fig. |U shows, for these relatively large values of a, BR(fi — > e + 7) 
and BR(r — > e + 7) can also be suppressed as a consequence of a partial cancellation between 
the contributions due to the complex R (r ^ 0) and the contribution present for R = 1 (see 
eqs. (jHEl)-(I2El))- For the chosen values mi and /3m, the effect of complex R on BR(r — > ^ + 7) 
is minor: the latter can change at most by a factor ~ 2 for r^5x 10 -2 . 

In the case of \a\, \b\, \c\ < 0.10 and tan/3 = 10; 50, illustrated in Fig. El r can take 
somewhat larger values. For tan/3 = 10, the contributions due to the matrix R in BR(/i — > 
e + 7) and BR(r — > e + 7) become dominant at m 1 ^ 0.01 eV; for BR(r — > fi + 7) the same 
is valid at m\ ^ 0.10 eV. Correspondingly, for m x ^ 0.01 eV (mi ^ 0.10 eV), BR(/i — > e + 7) 
and BR(r — > e+7) (BR(r — > /i+7)) can be enhanced significantly independently of the value 



7 We should note that some of the predicted values of the branching ratio BR(/i — > ej), presented in 
Figs. |3 and |5l arc already ruled out by the existing experimental limits on BR(/i — > e^). This illustrates, in 
particular, the type of constraints one might obtain, e.g. on the Majorana phase a and/or mi for a given 
specific set of values of SUSY breaking parameters. 
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of the Majorana phase a (Fig. EJ). For relatively large values of a, e.g. a ~ [7r/2, 7r], there is 
again the possibility of a reduction (or suppression) of the ij — > lj + 7 decay rates of interest 
with respect to those predicted for R = 1: the presence of the phase a can lead to partial 
cancellations between the three new terms proportional to \Jra\ra%, y/mim^ and ^Jrn^mi in 
eqs. Similar conclusions are valid in the case of tan/3 = 50. For 7r/2 ^ a < n 

in this case, the enhancement of BR(yU — > e + 7) and BR(r — > e + 7) at mi ^ 0.02 eV 
due to the RG effects is so strong that the contributions due to the complex R can lead to 
further increase of BR(ii — > e + 7) and BR(r — > e + 7) by not more than approximately 
an order of magnitude. For sufficiently large a, e.g. a ~ [tt/2, 7r], these contributions can 
also compensate partially the RG effect of enhancement leading to significantly smaller 
BR(/i — > e + 7) and BR(r — * e + 7) than in the case of R = 1 (Fig. EJ. Qualitatively similar 
results are valid also for BR(r — > // + 7). 

In the preceding discussion we have focused on the role of the Majorana phase a, ignoring 
the possible effects of the second Majorana phase in the PMNS matrix (3m- Figures |U and El 
correspond to (3m = 0. It follows from eqs. (|35 |) - (J35j) that if the phase (3m 7^ is sufficiently 
large, it can have an effect on the contributions to BR(/j — > lj +7) due to the complex matrix 
R, eq. (jnnj), similar to that of the phase a discussed above. 

6 Conclusions 

In the present article we have investigate the effects of the Majorana CP-violation 
phases in the PMNS matrix and of the RG evolution of the light neutrino mixing parameters 
- masses and PMNS mixing angles, on the predictions for the rates of lepton flavour violating 
(LFV) decays, fi — > e + 7, r — > e + 7 and r — > /i + 7, in a class of SUSY theories with see- 
saw mechanism of neutrino mass generation. We have considered minimal supersymmetric 
extensions of the Standard Theory with heavy RH Majorana neutrinos (MSSMRN), in which 
soft SUSY breaking by flavour-universal terms in the Lagrangian occurs at a scale Mqut 
above the RH Majorana neutrino mass scale Mr. In this class of theories, the charged lepton 
radiative decays yU^e + 7, r— ► /x + 7, etc. and other LFV processes are predicted to take 
place with rates that typically are within the reach of present and planned experiments. The 
predictions for the /i — > e + 7, r^ii + 7, etc. decay rates are known to depend critically 
on the matrix of neutrino Yukawa couplings at the scale M^j, Y^(M^). The matrix Y^(M^) 
can be expressed in terms of the light and heavy RH Majorana neutrino masses, the PMNS 
mixing matrix U - all at the scale Mr, and an orthogonal matrix R. The same matrix 
~Y u (Mr) is one of the basic ingredients of the see-saw mechanism and of the leptogenesis 
which can take place only if the orthogonal matrix R is complex. 

The neutrino mixing parameters entering into the expression for Y u (Mr) should be 
evaluated at the scale Mr. These include the solar and atmospheric neutrino mixing angles 
612 and 823, the CHOOZ angle #13, the two Majorana CP-violating phases in the PMNS 
matrix a and (3m ( see ec f ®)> an d the light neutrino masses m 1; m 2 = yjm\ + Am^, 
and m 3 = \Jm\ + Am^, where Am^ and Am^ are the neutrino mass squared differences 
driving the solar and atmospheric neutrino oscillations. In order to find the changes of 
the mixing angles 6^- and masses mj, due to the RG evolution from the scale Mz ~ 100 
GeV, where they are measured, to Mr, we have solved numerically the 1-loop RG equations 
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in MSSMRN for the see-saw Majorana mass matrix of the left-handed flavour neutrinos, 
eq. (|12j). The light neutrino mass spectrum was assumed to be with normal hierarchy, 
rrii < m 2 < m 3, and the lightest neutrino mass m x was allowed to take values from to 0.30 
eV. The current best fit values of the known neutrino mixing angles and AmL 31 ) at Mz 
(see eqs. © and ©) were used as initial conditions for the solution of the RG equations. 
We have presented results in the case of the CHOOZ mixing angle 8 13 and the Majorana 
phase /3m set to at Mz- In these and the subsequent calculations the scale Mr was fixed at 
Mr = 2 x 10 13 GeV. Our results concerning the RG evolution of the neutrino mixing angles 
and neutrino mass differences (m 2 (3) — mi) and Am^j) are discussed in Section HI and are 
illustrated in Figs. ^ and |2l where the indicated quantities are shown at the scale Mr as 
functions of m\. 

The values of the neutrino mixing parameters at the scale Mr were used as input in 
the calculation of Y U (M R ) and correspondingly of the rates of the LFV decays li — > lj + 7. 
The latter were performed for a number of typical sets of values of the soft SUSY breaking 
parameters (mo, m^, Aq) in the few 100 GeV range and for values of tan f3 from the interval 
(10-50). For simplicity, the analyses were performed for the case of quasi-degenerate in mass 
heavy RH Majorana neutrinos, Mi = M 2 = M 3 = Mr. The rates of the LFV decays 
/i — > e + 7, T^e + 7 and r — > /i + 7 have been calculated both in the cases of R = 1 and 
of complex matrix R. 

Our results show that if R = 1, taking into account the effects of the RG running of 
the neutrino mixing parameters in the case of relatively large tan/3 and/or mi, can change 
significantly the predictions for the li — > lj + 7 decay branching ratios BR(/i — > lj + 7). 
The magnitude of this change depends strongly on the value of the Majorana CP-violation 
phase a. More specifically, we find that for 0.05 ^ m x ^ 0.30 eV and 30 ^ tan/3 ^ 50, the 
predicted branching ratios BR(/i — * e+7) and BR(r — * e + 7), can be enhanced by the effects 
of the RG running of 9ij and (m.2(3) — mi) by 1 to 3 orders of magnitude if n/4 ^ a ^ tt, 
while BR(r — > fi + 7) can be enhanced by up to a factor of 10. For a = 0, the combined 
effect of the RG running of 612, #23 #13 and rrij on BR(/j — > lj +7) was found to be relatively 
small (essentially negligible) even for rrii ~ (0.2 — 0.3) eV and tan/3 ~ 50. 

We have considered further the case of non-trivial complex matrix R having the form 
given by eq. ()35|) . In this case R can be parametrised by three real parameters a, b, c. The 
form of R we have employed is the most general relevant one for the study of the LFV 
decays — > lj + 7 in the case of QD heavy RH neutrinos. As it has been shown in [28], 
successful leptogenesis requires abc 7^ and 10~ 5 ^ \abc\ 1. We find that in the case 
of tan/3 ~ 10 and mi ~ 0.06 eV, i.e., when the RG effects are negligible, the branching 
ratios BR(/x — > e + 7) and BR(r — > e + 7) can be significantly enhanced for any value of the 
Majorana phase a by the effect of complex R [28], provided r = \J a 2 + b 2 + c 2 ^ 5 x 10 -3 
(Fig. HJ) . The corresponding enhancement factor can be as large as ~ 100 for r = 5 x 10~ 2 . 
For relatively large values of a, e.g. a ~ [7r/2,7r], however, and depending on the signs and 
values of the parameters a, b, c, BR(/i — * e + 7) and BR(r — > e + 7) can also be suppressed as 
a consequence of a partial cancellation between the different contributions in the amplitudes 
of these two processes. The effect of complex R on BR(r — ► /i + 7) is minor for r ^ 5 x 10 -2 
and m x < 0.10 eV. In the case of a,b,c satisfying |a|, \b\, \c\ < 0.10, the contributions due 
to R in BR(yU — > e + 7) and BR(r — >• e + 7) become dominant at mi ^ 0.01 eV; for 
BR(r — > ji + 7) they become dominant at mi ^ 0.10 eV. Similar results are valid in the 



15 



case of tan/3 ~ 50. For n/2 ^ a < it in this case, the enhancement of BR(/i — > e + 7) and 
BR(r — > e + 7) at m x ^ 0.02 eV due to the RG effects is so strong that the contributions of 
the complex R can lead to further increase of BR(/i — > e + 7) and BR(r — > e + 7) by not 
more than approximately an order of magnitude. For sufficiently large a, e.g. a ~ [7r/2,7r], 
these contributions can also compensate partially the RG effect of enhancement leading to 
significantly smaller BR(/i — > e + 7) and BR(r — > e + 7) than in the case of R = 1 (Fig. EJ). 
Qualitatively similar results are valid also for BR(r — > yU + 7). If the Majorana phase (3m is 
sufficiently large, it can have an effect on the contributions to BR(/j — > Lj + 7) due to the 
complex matrix R, similar to that of the phase a discussed above. 

The results of the present study show that taking into account the effects of the Majorana 
CP- violation phases in the PMNS matrix and of the RG evolution of the light neutrino mixing 
parameters - masses and PMNS mixing angles, can have a remarkably large impact on the 
predicted rates of the LFV decays /1 — > e + 7 and r — > e + 7, etc. in MSSM with right-handed 
neutrinos and see-saw mechanism of neutrino mass generation. The impact of these effects 
on the r — > fi + 7 decay rate is typically smaller and possible only for mi ^ 0.10 eV. 
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APPENDIX A 

In the present Appendix we derive approximate analytic expressions for some of the 
neutrino mixing parameters at the scale Mr. The renormalisation group (RG) equation for 
the neutrino mass matrix m u in MSSMRN has the form (see, e.g., [37]) 



d 



din fx R 



m. 



8vr 2 



4tt Ua 2 + jjo^ + tr(3y^)| m v + l - {(Y ] E Y E ) T m v + m u (Y^Y E )} 



(39) 



where fiR is the RG energy scale, ctj = c/|/(47r), where #3, g 2 and g\ are the SU(3) C , SU{2)l 
and U(1)y gauge coupling constants, respectively, and Yjj is the matrix of Yukawa couplings 
of up-quarks. This equation admits the approximate solution [34] 



m v [M R ) ~ J (M R ) 



where J(M R ) and e r are given by 



m v {M z 



(40) 



1 + e 7 



J (M R ) =exp 
e T =l- 



2 

M z 



In M R 



n JlnMz 

(l/87r 2 )(l+tan 2 /3)(m T /u) 2 



tr(3Y^)Yu(fi)) - 4tt ( 3a 2 (/x) + ^(/i) ) ) d(lnji) 



3 
— < 

5 



M, 



(41) 
(42) 
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m T being the r mass. Note that e r takes positive values and is a monotonically increasing 
function of tan/5. The PMNS matrix at M R , U(M R ), is obtained by diagonalising m u (M R ): 

U(M R fm u (M R )U(M R ) = diag(m 1 (M R ),m 2 (M R ),m 3 (M R )) . (43) 

Using eqs. (jU) and (f4T)j) - (fHIj) it is possible to derive the following approximate expressions 
for the solar neutrino mixing angle 9\ 2 and Am^ at M R . 

|sine 12 (M z )cos^. + cos^ 12 (Af z )sin^.e- iQ / 2 | 
tanMM*)^ | co ^ 12(Mz )cos^-sin^ 12 (M z )sin^e-W 2 | ' W 

and 

A^ l( M H ) s J\M R ) Am "< Mz > + 4f - m ' (Mz) 3i " 2 ^' il -fe) , (45) 

cos 26^ 

where the angle 9 X is determined by the equation 

4e r m 2 (M z ) sin 2 6 23 (M Z ) sm26 12 (M z ) cos ^ 

tciii 2aj — — ( 4o ) 

Am 2 21 (M z ) +4e T ml(M z ) sin 2 9 23 (M Z ) cos 29 12 (M Z ) ' V ; 

As can be shown, eqs. ff43j) and (j4T)j) imply cos 2^ > 0, sin 2^ < 0, provided we keep the 
ordering mi(M R ) < m 2 (M R ) of the masses mi and m 2 at M R . It follows, in particular, from 
eq. (jSj) that 

sin 2 29 12 (M R ) < sin 2 29 12 (M Z ) . (47) 
For a(M z ) = and a(M z ) = it, for instance, we have 



sin 2 2(MM z ) + ^.) , « = 
sin 2 29 12 (M z ) , a = vr. 



sin 2 2^ 12 (M i? ) = <!^_ 2o y;^r " ' (48) 



The mixing angle #13 and the Dirac CP-violation phase 6 are also affected by the RG 
evolution from M z to M R . If, e.g. we set s\ 3 = at Mz, we get for S13 and <5 at M^: 

J 2 (A.f fl )e r m 1 (M 2: )m 3 (M z )sm2e 12 (M 2; )sin29 2 3(Af i; ) . a(M z ) 
W " ml(Afe) - m ; ( Af fi ) S1 " — ' (49) 



and 



where 



5(M ii )~ei + e2-| + ^^-^(M^), fora(M z )^0, (50) 



6 = arg(c 12 (M z ) Cl . - Sl2 (M z ) Sx e-^ M ^ 2 ) , 

6 = arg( Sl2 (M z )c x + c 12 (M z ) S:c e^ M ^/ 2 ) . (51) 

The angle #13 (Mr) takes its largest value for a(M z ) = it and in this case S(M R ) = —/3m(M z ). 
In contrast, if a = 0, the correction to sin 9 13 due to the RG running is negligible. 
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Figure 1: The dependence of tan 2 6* sol = tan 2 # 12 , tan 2 atm = tan 2 623 and of | V13 j = sin 6*13, 
evaluated at Mr, on the lightest neutrino mass ni\ for tan/3 = 10 and tan/3 = 50. The 
following values neutrino mixing parameters at Mz were used as input in the calculation 
of the RG running effects: Am| = 8.0 x 10~ 5 eV 2 , Am 2 A = 2.2 x 10~ 3 eV 2 , tan 2 6 Q = 0.4, 
tan 2 9a = 1, and sin^ 13 = 0.0. The neutrino mass spectrum at M z is assumed to be with 
normal hierarchy, m\{Mz) < ^{Mz) < m^Mz)- The RH neutrino mass scale Mr is fixed 
as M R = 2 x 10 13 GeV. 
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Figure 2: The light neutrino mass differences at the scale Mr as a functions of m\ for 
a = 0;tt and tan/? = 10; 50: (a) Am 2 A (M R ) = m\ (Mr) - m\{M R ), (b) Am%(M R ) = 
ml(M R ) - ml(M R ), (c) m 2 {M R ) - mi(M R ), and (d) m 3 (M R ) - mi(M R ). The neutrino 
mixing parameters at Mz and the mass M R used are the same as in Fig. 
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Figure 3: The branching ratios of the LFV decays yU — > e + 7, + 7 and r — > e + 7 

versus mi in the case of R = 1 and for a = 0; 7r/2; 7T. The results shown correspond to 
the following two sets of SUSY parameters: (A) tan/3 = 10, m = 100 GeV, m 1 / 2 = 250 
GeV, A = -100 GeV, and (B) tan/5 = 50, m = 400 GeV, m l/2 = 300 GeV, A = 0. The 
neutrino mixing parameters are the same as in Fig. ^ 
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Figure 4: The branching ratios BR(/i — > e + 7), BR(r — > e + 7), and BR(r — ► + 7) as 
functions of r = + & 2 + c 2 for fixed mi = 0.06 eV, a = 0; 7r/2; 7t and = 0. The 
values of the other neutrino mixing parameters at Mz are the same as in Fig. ^ The SUSY 
parameters are taken as tan/3 = 10, muz = 250 GeV, m = 100 GeV and A = —100 GeV. 
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Figure 5: The same as in Fig.|3]but for complex matrix 1, eq. (jSHJ)- The three constants 
a,b,c parametrising the matrix R are assumed to lie in the interval —0.1 < a,b,c < 0.1. For 
comparison results for real R = 1 are also shown (see text for further details). 
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